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Electron Confinement, Orbital Ordering, and Orbital Moments 
in (f-d^ Oxide Heterostructures 

Victor Pardo^'^'B and Warren E. Pickett^'Q 
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The (SrTi03)m/(SrV03)n d^ — d^ multilayer system is studied with first principles methods 
through the observed insulator-to-metal transition with increasing thickness of the SrVOs layer. 
When correlation effects with reasonable magnitude are included, crystal field splittings from the 
structural relaxations together with spin-orbit coupling (SOC) determines the behavior of the elec- 
tronic and magnetic structures. These confined slabs of SrVOs prefer Qorb={7r,7r) orbital ordering 
of iz =0 and iz = —1 {jz = — ^) orbitals within the plane, accompanied by Qspin={Ofi) spin 
order (ferromagnetic alignment). The result is a SOC-driven ferromagnetic Mott insulator. The 
orbital moment of 0.75 /is strongly compensates the spin moment on the iz = —I sublattice. The 
insulator-metal transition for n = 1 ^ 5 (occurring between n=4 and n=5) is reproduced. Unlike 
in the isoelectronic d*^ — d^ Ti02/V02 (rutile structure) system and in spite of some similarities in 
orbital ordering, no semi-Dirac point [Phys. Rev. Lett. 102, 166803 (2009)] is encountered, but 
the insulator-to-metal transition occurs through a different type of unusual phase. For n=5 this 
system is very near (or at) a unique semimetallic state in which the Fermi energy is topologically 
determined and the Fermi surface consists of identical electron and hole Fermi circles centered at 
A:=0. The dispersion consists of what can be regarded as a continuum of radially-directed Dirac 
points, forming a "Dirac circle". 



I. BACKGROUND 

Oxide heterostructures with a polar discontinuity 
across interfaces have attracted a great deal of atten- 
tion recently, due to the unusual electronic behavior 
that can arise. ^^^ It is now becoming evident that 
heterostructures with non-polar interfaces can also 
lead to unanticipated behavior, including low en- 
ergy dispersion that is distinct from any previously 
known system. The specific example is the dP/d^ 
interface system Ti02/V02 that displays a point 
Fermi surface, from which semi-Dirac dispersion^ 
emerges. Semi-Dirac dispersion is characterized by 
conventional, massive dispersion along one direction 
in the two-dimensional plane but massless dispersion 
in the perpendicular direction. 

Although reminiscent of graphene, the semi- 
Dirac system displays its own distinctive low en- 
ergy properties.^ The behavior is actually an elec- 
tron confinement phenomenon assisted by a partic- 
ular orbital ordering, and these nanostructures dis- 
play a peculiar metal-insulator transition^ as the 
thickness of the d^ oxide is increased. In these 
rutile structured oxides, the metal- insulator transi- 
tion takes place through an intermediate semi-Dirac 
point when the thickness is approximately 1 nm, 
where the system is neither insulating nor conduct- 
ing and the Fermi surface is point-like. 

Transition metal oxide perovskites with 3d^ con- 
figuration are known to be on the borderline between 
metallic and insulating, depending on the relative 
sizes of several electronic energy scales, including 
the ratio between the on-site Coulomb repulsion U 



and the bandwidth W of the d electrons, and the 
competition between magnetic energies and Jahn- 
Teller splittings. For small U/W, the material will 
be metallic, like the correlated metal^ SrVOs. For a 
large value of the ratio U/W, the system will present 
a more localized behavior, and a Mott insulator will 
result, as in LaTiOa or YTiOa.^ Other energy scales 
may also affect, or even determine, the delicate bal- 
ance. 

Multilayers of SrVOs (SVO) and SrTiOa (STO) 
have been grown on STO substrates by Kim et at. 
and have displayed a transition^ from the typical in- 
sulating behavior of STO to the metallic behavior 
of SVO as the number of layers of each constituent 
is increased. Superlattices formed by films with 2 
layers of STO and 6 layers of SVO [we will denote 
(SrTi03)m/(SrV03)n (001) oriented multilayers as 
m/n] already show metallic behavior with a nearly 
flat resistivity curve with magnitude close to that 
of SVO. However, from the behavior of the resistiv- 
ity, it was observed that the 2/3, 2/4 and 2/6 films 
(increasing SVO thickness) show an insulator-metal 
transition at temperatures ranging from about 100 
K for the 2/6 film to the approximately 230 K of 
the 2/4 and 2/3 films. It is to be noted also that 
only 2 layers of STO are not enough to isolate the 
SVO slabs, as can be seen from the resistivity data 
for the 2/3 and 3/3 systems. In fact, the 2/3 system 
is already semi- metallic, but four layers of STO are 
enough to render interactions between neighboring 
SVO slabs negligible. Five SVO layers are needed to 
obtain a metallic state, as we show below from our 
calculations. 



In this paper we extend our investigations of dP-d^ 
nanostructures by studying this STO/SVO system, 
focusing on the differences that the crystal structure, 
with its specific crystal field splittings, can cause. 
We choose the most commonly studied structure, 
perovskite, with the previously studied nanostruc- 
tures with rutile crystal structure. We will compare 
different orbital orderings and magnetic arrange- 
ments, and also study how the insulator to metal 
transition occurs with increasing thickness of the d^ 
material (SrVOa). 



II. COMPUTATIONAL METHODS 

Our electronic structure calculations were per- 
formed within density functional theory^ using the 
all-electron, full potential code wien2k^^ based 
on the augmented plane wave plus local orbital 
(APW+lo) basis set^^. The exchange-correlation 
potential utilized for the structure optimizations was 
the generalized gradient approximation (GGA) in 
the Perdew-Burke-Ernzerhof (PBE)^^ To deal with 
strong correlation effects that are evident in SrVOs 
we apply the LSDA+U scheme^^^^^ including an on- 
site U and J (on-site Coulomb repulsion and ex- 
change strengths) for the Ti and V 3d states. The 
values U= 4.5 eV, J= 0.7 eV have been used for V to 
deal properly with correlations in this multilayered 
structure; these values are comparable to those used 
in literature for d^ systemsi^"— . Our calculations 
show that a larger U, above 5.0 eV, give an incorrect 
insulating behavior of bulk SrVOs in cubic struc- 
ture, hence overestimating electron-electron interac- 
tions. Since Ti d states never have any significant 
occupation, including U or not on Ti 3d orbitals has 
negligible consequence. Spin-orbit coupling (SOC) 
effects on the valence and conduction states, which 
is discussed in the last section, have been introduced 
in a second variational procedure built on the scalar 
relativistic approximation. 



III. TREATING CORRELATION EFFECTS 
IN BULK SRVO3 

We have first performed calculations on bulk 
SrVOs (SVO) to establish the thick SrVOs limit 
of these nanostructures. Experimentally, SrVOs is 
a ferromagnetic (FM) metal^^ crystallizing in a cu- 
bic perovskite structure. ^^ No distortion from cubic 
structure has been observed experimentally, which is 
consistent to its metallic character; a Mott insulat- 
ing d^ system would be expected to distort due to or- 
bital ordering. Our calculations show that the most 
stable structure based on GGA exchange-correlation 



(not including on-site Coulomb repulsion effects) is 
cubic with no distortions. However, a non-magnetic 
solution is obtained as a ground state within GGA. 

When correlations are introduced by means of 
the LDA+U method, the most stable structure is 
slightly distorted. This broken symmetry arises be- 
cause the LDA+U method tends to promote integer 
occupations of one of the t2g orbitals (which becomes 
lower in energy, and preferentially occupied). 

Calculations were carried out with a tetragonal 
distortion in bulk SrVOs and with various values of 
the on-site Coulomb repulsion U. LDA+U calcula- 
tions predict two possible orbital configurations: a 
FM configuration with all the t2g bands equally oc- 
cupied and an AF solution with a substantial occu- 
pation of the dxy orbital. As U increases above 5 eV, 
the AF solution is more stable, even for undistorted 
cubic perovskite, leading to an incorrect AF insulat- 
ing state. In all ranges of U studied (3 to 7 eV), a 
tetragonal distortion is more stable, whereas a sim- 
ple GGA calculation leads to an undistorted cubic 
solution, in agreement with experiment. Even in the 
case of a FM solution, when a U is introduced in the 
calculations, a tetragonal distortion is obtained as 
a ground state, in disagreement with experimental 
observations. Using the correct structure obtained 
from both experiment and GGA calculations (which 
agree), the LDA+U method with values of U smaller 
than 5 eV leads to the correct FM metallic state as 
ground state. 

Within GGA (U = 0) SVO is a cubic metal, with 
all the t2g orbitals are equally occupied. Within 
LDA+U, a tetragonal distortion causes the prefer- 
ential occupation of the dxy orbital, an insulating 
state and also an antiferromagnetic (AFM) ordering 
to be stabilized, whereas a cubic FM metallic state is 
observed experimentally. This is the case also when 
LDA+U is applied for the structure relaxed with 
GGA. For this reason, throughout the paper, we will 
use for structure optimizations of these multilayers 
including such a moderately correlated compound, 
the GGA-PBE functional J^ With the structure thus 
determined, we use the LDA+U method for the cal- 
culations of the electronic and magnetic structure, 
and energy differences. In bulk SVO this procedure 
properly results in a FM half metal with 1 /i^/V 
atom. 



IV. GENERAL CONSIDERATIONS 

SrVOs has a lattice parameter of 3.84 A,^^ while 
the lattice parameter of SrTiOs (STO) is 3.90 A.^ 
Since most of these superlattices are grown on a 
SrTiOs substrate, for closest comparison with exper- 
imental data we fix the a lattice parameter to be the 



STO lattice parameter. However, we use the GGA- 
PBE exchange-correlation potential to optimize the 
superlattice c parameter and also the internal co- 
ordinates of all the atoms by a force minimization 
together with a total energy minimization. 

Since the IF between SVO and STO has no polar 
discontinuity, the distortions introduced at the inter- 
face between the two oxides by the nanostructuring 
can be understood as first, strain due to the lattice 
parameter mismatch (1.5%), and secondly, to charge 
imbalance within the V t2g orbit als. However, no 
ionic charge compensation effects of the sort that 
are so interesting in LaAlOs/SrTiOa nanostructures 
are present i*^"— 

We compare results for m/n multilayers with m= 
4 layers of STO (about 1.6 nm thickness) sandwich- 
ing an SVO layer with variable thickness from 1 to 5 
layers of SVO (0.4 to 2.0 nm), because we find 4 lay- 
ers of STO sufficient to isolate the SVO slabs to give 
two-dimensional behavior (negligible k^-dispersion 
in the band structure). We analyze the evolution 
of the electronic structure for increasing thickness of 
SVO slabs (n from 1 to 5). 



A. d^ V ion in an octahedral environment 




FIG. 1: (Color online) Spin density isosurface of the n—2 
system at 0.8 e/A^. Shown is a top view of the x — y 
plane in the planar AF Mott insulator state, which is 
the most stable only at larger values of U. The da:;^ or- 
bital is occupied on each V ion. 
different spin directions. 



Different colors indicate 



One common feature within the SVO sublayer for 



all thicknesses we have studied is that the V octa- 
hedral environment, especially at the interface (IF), 
will be tetrahedrally distorted from its cubic sym- 
metry in the bulk. It was noted above that such a 
distortion, treated fully (i.e. including structural re- 
laxation) within the LDA+U method for reasonable 
values of U, produces an (incorrect) AF insulator 
state for bulk SVO. In the case of the multilayers, 
our relaxation of the c lattice parameter and atomic 
positions will lead to values of the interplane V-V 
distance somewhat smaller than the in-plane value. 
The simplest scenario would be that the multilayer 
structuring of SVO on an STO substrate will pro- 
duce a tetragonal distortion of the oxygen octahedra 
around the V cations. This c-axis contraction of the 
octahedra leads to a preferred occupation of the dxy 
orbital. If dxy orbitals are occupied on all atoms in 
a layer, a small enough bandwidth or strong enough 
intra-atomic interaction U will give AF ordering by 
superexchange. This effect is observed experimen- 
tally in similar multilayer systems with the d^ com- 
pound being a Mott insulator like LaTiOs-^^ 

A more realistic picture must allow for the likeli- 
hood that the degenerate dxz-,dyz orbitals will have 
a narrower bandwidth in the x — y plane, and cor- 
relation effects may favor occupation of some com- 
bination of these orbitals. Even though small, spin- 
orbit coupling breaks this degeneracy and specifies a 
favored combination. The actual occupation will de- 
pend on several factors. First of all is the tetragonal 
crystal field splitting of the t2g orbitals. In all cases 
we study, the xy on-site energy is lower due to in- 
plane stretching imposed by the STO substrate. Sec- 
ond, the X — y plane bandwidths of dxy and dxz^ dyz 
bands are very different, and will depend on orbital 
order. Thirdly, the orbital order is closely tied to the 
magnetic configuration of the system. The simplest 
possibility (above), with all the electrons in an dxy 
orbital, leads to AF order (this orbital pattern can 
be seen in Fig. [1]), whereas we find that occupying 
dxzidyz orbitals favors FM order. 



B. Structural relaxation: strain effects 

The structural distortion as n increases can be de- 
scribed in terms of lattice strain along the c-axis. Re- 
laxing the c-axis value (and also the internal atomic 
positions) yields the results provided in Table [H For 
understanding the structural distortions, we can de- 
fine four different distances along the c-axis (within 
the plane they are constrained by the STO lattice 
parameter 3.905 A): the V-V distance, the V-Ti dis- 
tance across the IF, the Ti-Ti distance, and finally 
the average c lattice parameter. We find that the 
Ti-Ti distance hardly changes, and also the V-Ti 



distance variation is minor. However, the V-V dis- 
tance grows towards a limiting value as the number 
of SVO layers increases. Both V-V and Ti-Ti in- 
terlayer distances presented in the table vary only 
by ± 0.01 A for the various layers within each n/m 
system. 



TABLE I: Interlayer distances (A) between metal 
cations, for the various n/m configurations under study. 
The important distances to describe the structure are Ti- 
Ti, Ti-V (across the interface), V-V, and the average c 
parameter. The Ti-Ti separations, and V-V separations, 
for each heterostructure are uniform to within 0.01 A. 





Ti-Ti 


V-Ti 


V-V 


^av 


4/1 


3.97 


3.92 


- 


3.95 


4/2 


3.95 


3.90 


3.79 


3.91 


4/3 


3.95 


3.91 


3.83 


3.91 


4/4 


3.98 


3.94 


3.86 


3.92 


4/5 


3.98 


3.93 


3.87 


3.92 



C. Band lineups; intralayer supercells 

Some general features should be established. Due 
to the alignment of the SrVOa Fermi level (or gap) 
within the SrTiOa gap between filled O 2p states 
and empty Ti Sd states, there is a 2.5 eV energy 
window which contains only the V d bands that are 
of interest for us. The oxygen 2p bands lie below the 
V d bands. We note especially that we have used 
laterally enlarged, c(2 x 2) superstructures to allow 
the possibility (or likelihood) of AF spin alignment 
as well as FM, and various orbital ordering patterns 
as well. The ground states we obtain and analyze are 
all lower in energy than any that would have been 
found in primitive p(l x 1) cells. 



D. Spin-orbit coupling 



the orbital moments are typically quenched by mix- 
ing with neighboring orbit als. Fairly recently sev- 
eral cases have come to light in which t2g orbital 
moments can be quenched^^ by structure-induced 
or spontaneous symmetry lowering, or even in a bd 
system an orbital moment can compensate a spin 
moment and prevent orbital ordering. ^^ The calcu- 
lations presented here provide an additional example 
of the importance of SOC in Zd systems. 

The strain in the slabs we study breaks the t2g 
symmetry, with dxy lying lower. Strong magne- 
tocrystalline anisotropy favors the magnetization 
along the c-axis (by 130 meV/V compared to an 
in-plane orientation), making this the quantization 
axis. Then one can expect competition between oc- 
cupation of the d- and dxy orbitals in a d^ ion, with 
kinetic energy (bandwidth and band placement) be- 
ing a determining factor. The energy gain, and thus 
the large magnetocrystalline anisotropy, is related 
(see following sections) to the formation of a large 
orbital moment 0.75 jib along the z-axis when the 
d- orbital is occupied. Since the spin moment is re- 
duced somewhat from its atomic value of 1 iib by 
hybridization, the net magnetic moment on such an 
ion can be quite small. 

In the next section we show that an alternating 
orbitally-ordered (AOO) state, with half the elec- 
trons in a dxy orbital and the other sublattice in a d± 
orbital, arises and leads to FM spin alignment. The 
AF solution becomes favored only at unphysically 
large values of U (above 5-6 eV). In the reasonable 
range of values of U, this AOO state is energetically 
favored for several SVO slab thicknesses. This AOO, 
Mott insulating FM state competes and overcomes 
the Mott-insulating AF state, which has all the elec- 
trons in xy orbitals. Unexpectedly, when SOC is 
included, the alternating orbital ordering produces 
an FM Mott insulating state at small SrVOa thick- 
nesses, below 2 nm, and for realistic values of U for 
this multilayered system. It is instructive to follow 
the behavior through the insulator-to-metal transi- 
tion with SVO thickness. 



SOC in the V atom usually produces minor effects, 
except for those which depend entirely on it (such as 
magnetocrystalline anisotropy). We find here how- 
ever, irrespective of the number of VO2 layers, that 
SOC in conjunction with correlation effects (the on- 
site repulsion U) completely alter the ground state 
that we obtain. 

It has long been known that the t2g subshell pro- 
vides a representation for ^=1 (not ^=2) orbital 
moments^ -^^ The rui states are d^ = dxy^ and d± = 
dxz ^idyz. SOC splits these, with d- {jz = — ^) ly- 
ing lower. This viewpoint is not common, because 



V. EVOLUTION OF THE ELECTRONIC 
STRUCTURE 

A. n= 1 confined layer of SrVOs 

When a single layer of SVO is confined by insulat- 
ing STO, the resulting strain lowers the dxy orbital 
energy, and it should be expected that a candidate 
ground state is dxy "orbitally ordered" AF insula- 
tor due to superexchange. Indeed this Neel state 
can be obtained in our calculations, for the moder- 
ate range of U that is relevant. The coupling arises 
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FIG. 2: (Color online) Majority spin LDA+U band structure of the n=l system without (left panel) and with 
spin-orbit coupling (right panel) included. A SOC-driven FM Mott insulator regime occurs at this very small SVO 
thickness: the introduction of SOC breaks the dxz^dyz symmetry, and LDA+U splits the d- (down) and d-\- states 
(up), opening the gap. The lower and upper Hubbard bands lie at -0.5 eV and 1.3 eV, respectively. 



through ddiT-type hopping between dxy orbitals in 
the plane (see Fig. [1]) as expected from Goodenough- 
Kanamori-Anderson (GKA) rules. ^^ However, the 
AGO phase with FM spin alignment was also ob- 
tained, and it is lower in energy by 4 meV/V. The 
two sublattices lead to distinct sets of bands, as eas- 
ily seen from the left panel of Fig. [2l The dxy 
band has the familiar square-lattice shape (distorted 
somewhat by 2nd neighbor coupling, and is 1.5 eV 
wide. The dxz — dyz bands are much narrower (0.25 
eV) because there is only ddir coupling. The centroid 
of the latter pair lies about 0.5 eV above that of the 
dxy band, providing the magnitude of the splitting 
of the t2g by strain. The system is metallic, with all 
three bands leading to Fermi surfaces. 

The picture is completely changed by SGC. Figure 
[2] compares the majority spin band structure of the 
FM AGG state, first without SGC, then with SGC 
included. SGC has no effect on the dxy {iz = 0) 
band. However, SGC breaks the symmetry of the 
iz = ^l doublet, and the narrow bandwidth (0.25 
eV) compared to the value of U results in a Mott 
insulating type of splitting of the d- and d-^ bands, 
by roughly ±U. The SGC-driven symmetry lowering 
is leveraged by the strong on-site interaction. The 
result is an AGG FM Mott insulator with a gap of 
0.6 eV. This d- orbital acquires a large orbital mo- 



ment of about 0.75 /i^, strongly compensating the 
spin moment. This motif of FM AGG V orbitals will 
recur for thicker SVG slabs. 



B. n= 2 SrVOs layers 

The bands near the gap for the two SVG layer slab 
are displayed in Fig. [3] without SGC included, to 
illustrate that for two layers there is already a band 
gap without SGC, produced by interlayer coupling 
of dxz ^ dyz states. However, SGC produces the same 
orbital ordering and intralayer FM alignment as for 
n=l, and the layers are also spin- aligned to give an 
overall FM AGG Mott insulating state. Comparing 
the total energy for the dxy AF state with the more 
stable FM state, the energy difference is found to 
be large, 76 meV/V. This energy difference includes 
the in-plane energy gain as for n=l and the larger 
interplane stabilization due to a large dda coupling 
between d_ orbitals. 



C. n= 3 SrVOs layers 

As mentioned earlier, structural relaxation was 
always carried out using only the GGA exchange- 
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FIG. 3: (Color online) LDA+U band structure (without 
SOC) of the 4/2 system (U = 4.5 eV). A Mott-insulator 
regime occurs at this thickness of two SVO monolayers. 



FIG. 4: (Color online) Band structure of the 4/3 system 
for U= 4.5 eV with spin-orbit coupling spin state. The 
minority spin state would be equivalent in this AFM 
ground state. A very small gap insulator is obtained, 
on the verge of a metal-insulator transition. A metal is 
actually obtained at small U (below 4 eV). 



correlation functional. We note that for this n=3 
case GGA calculations give an AFM interplane cou- 
pHng (up/down/up) of FM layers. We carried out 
the structural relaxations in this magnetic structure; 
however, the type of magnetic order is not expected 
to affect the relaxation appreciably. 

Returning to the LDA+U calculations to evaluate 
spin and orbital order, the ground state is FM overall 
(up/up/up) with AGO order within the plane, and 
a metallic band structure. As for n=l and n=2, in- 
plane AFM ordering with all dxy orbitals occupied 
can be obtained, but is energetically higher than the 
FM AGO state. In this case, we can compare the 
total energies of two configurations with the same in- 
plane orbital ordering, but different inter layer AGO 
alignment, AAA (like orbitals aligned along the c- 
axis) or ABA. The AAA configuration is more stable 
by 6 meV/V, giving an idea of the strength of the 
interlayer coupling. 

SGC again results in the same intralayer AGG 
FM state, however in this case there is an extremely 
small gap at T, visible in Fig. HJ It differs from 
the n=2 system (and is similar to the n=l system) 
in having a gap only when SGC in included (break- 
ing the d± degeneracy that otherwise leaves half (or 
partially) filled bands). 



D. n=4 SrVOs layers 

The band structure of the n=4 system, which has 
the same AGG FM Mott insulating ground state, 
is shown neglecting SGC in Fig. [5] to show that 
SGC is not necessary to produce a gap (as for n=2), 
although it does change the band structure. The 
four nearly degenerate dxy bands with their sim- 
ple square-lattice shape lie in the -1.6 to -0.1 eV 
range, interlayer coupling is very small for these dxy 
orbitals. The four band pairs for the dxz^dyz or- 
bitals correspond to linear combinations of similarly 
shaped bands for each of the four layers. The band 
dispersion decreases from 0.7 eV for the lowest band 
to zero for the highest (unoccupied) bands. By using 
a simple tight-binding model with interplane cou- 
pling, the band ordering and hence the occupation 
of the bands is readily reproduced. 



E. n=5 SrVOs layers 

For n=5 several orbital orderings can be obtained 
and their energies compared. Gf the solutions we 
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FIG. 5: (Color online) Majority spin band structure of 
the n=4 system (U= 4.5 eV) of this FM half- metallic, 
orbit ally ordered ground state. The four characteristic, 
nearly degenerate dxy bands lie in the -1.6 to -0.2 eV 
range. There is ferro-orbital ordering between planes. 
The inter-layer coupling splits the d- bands by 0.2-0.3 
eV. 







FIG. 6: (Color online) Spin density isosurface of the n=4 
system (U=4.5 eV) of the majority spin electrons of this 
FM insulating solution. The in-plane alternating orbital 
ordering dxy and d- orbitals is apparent, and also the 
interlayer orbital configuration. The stabilization of the 
ferro-orbital occupation along the z-axis is stabilized by 
the interaction between V dxz±iyz orbitals in this FM 
ground state. 



studied, again the AOO FM layers, spin aligned be- 
tween layers (globally FM) is again the ground state. 
Tetragonal strain and SOC result in the same AOO 
arrangement and corresponding identifiable bands, 
which are shown in Fig. [71 The five dxy bands have 
a total splitting no more than 0.4 eV, and are com- 
pletely filled. The five d- bands, each with small 
bandwidth, cover a range of 1.4 eV. The lower two 
(-1.4 eV to -0.9 eV) mix with the O 2p bands, and 
are not as obvious as the upper three bands. 

The distinction for n=5 is that the uppermost d- 
band is not fully occupied, but overlaps by 0.2 eV 
a conduction band that dips below the d- at the 
r point and leads to a semimetal. The insulator- 
to-metal transition has occurred between n=4 and 
n=5, very similar to the transition observed by Kim 
et al^ in their SVO/STO multilayers. Unlike the 
V02/Ti02 system where the transition proceeds 
through a point Fermi surface, semi-Dirac phased"— 
this transition appears to occur in the classic fash- 
ion of band overlap. The strong effects of tetrag- 
onal strain (symmetry breaking) and SOC (further 
symmetry breaking) are finally overcome by the in- 
creasing delocalization across the SVO slab as the 



quantum confinement effects are eroded. 



VI. TRANSITION THROUGH A 
SINGULAR FERMI SURFACE 

A very unconventional insulator-to-metal transi- 
tion nearly occurs at the n=5 thickness, and might 
actually occur for somewhat different value of intra- 
atomic repulsion U or different value of strain. If 
some small change lifted the conduction band at the 
M point in Fig. [7] above the Fermi level (it over- 
laps only slightly as it is), and if the overlapping va- 
lence and conduction bands at F are in the (small k) 
quadratic limit so their constant energy surfaces are 
circles, then the equality of electron and hole density 
leads to coinciding electron and hole Fermi surfaces. 
The "Fermi surface" is actually two identical elec- 
tron and hole Fermi lines, the "Fermi surface" has 
a boundary but no area. The bands near the Fermi 
energy, given in the simplest model by 



e/e = ±v{\k\ - kp), 



(1) 



are presented in Fig. [H where these linearly dis- 
persing bands are shown. With doping, the electron 
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FIG. 7: (Color online) Band structure of the 4/5 system 
(U=4.5 eV) of the majority spin of this FM half-metallic 
state. There is orbital ordering in-plane, and ferro- 
orbital ordering between layers, as for thinner SrVOs 
slabs. However, with the added bands, the interlayer 
coupling is no longer large enough to open a gap, leav- 
ing a semimetallic state with band overlap at F. 




FIG. 8: (Color online) Two (locally) linearly dispersive 
bands in the form efc= ± v||/c| — /cf| come together form- 
ing a continuous circle of Dirac points, i.e. a "Dirac cir- 
cle." Only one quadrant of the zone is displayed. This 
model represents the two bands closest to the Fermi level 
for the n= 5 system displayed in Fig. [71 



is even and odd symmetry under z-reflection, dxy 
states for example being even, and dxz^dyz states 
being odd. This is the case in which this Dirac- 
continuum of points may arise. 



and hole surfaces (lines, since this is two dimensions) 
separate, leaving an annulus that contains electrons 
if electron-doped, or holes if it is hole-doped. It is 
thus easy to understand how the annulus vanished 
as the doping level vanishes. The situation is in fact 
a continuum of radial Dirac points. Conversely, it 
presents the limit of a semi-Dirac point^^^ when the 
effective mass diverges. 

From Fig. [71 it is clear that (with the conduction 
band at M out of the picture) the positioning of the 
Fermi level at the band crossing point is topologically 
determined: only for precisely that Fermi level are 
an integral number of bands occupied, which is ex- 
actly what is required to occupy the V 3d electrons. 
The lowering of the symmetry of the eigenstates off 
of the symmetry directions (corresponding to some 
subgroup of the group of the symmetric k-point) may 
lead to coupling of the bands and the opening of a 
gap away from the X or M points (or both). In this 
case, one is left with four (respectively, eight) semi- 
Dirac points along symmetry lines. 

However, it may occur that the bands have dif- 
ferent symmetry throughout the Brillouin zone, in 
which case no gap opens. The simplest example 



VII. SUMMARY 



We have analyzed the transition from insulator 
to metal in the electronic structure of multilayers 
(SrTi03)4/(SrV03)n, with n varying from 1 to 5. 
The transition is observed to occur between n=4 and 
n=5. The origin of the observed changes is surpris- 
ingly intricate, with tetragonal strain and spin-orbit 
coupling (each with an associated symmetry break- 
ing) leveraging strong interaction effects modeled by 
the LDA+U approach. The effects of quantum con- 
finement finally determine the conduction charac- 
ter. Insulating behavior with a peculiar alternating 
orbital ordering within each V layer and FM mag- 
netic order results from a ferromagnetic Mott insu- 
lating state for n=4 or less. Ferromagnetic Mott 
insulators are rare, and these results indicate how 
this kind may be achieved (even designed) in ox- 
ide nanostructures. The FM insulator-to-FM metal 
transition finally results from band overlap as quan- 
tum confinement effects decrease. This system is 
very close to, if not at, an unusual semimetal state 
for n=5 in which the Fermi surface is topologically 
determined and consists of two degenerate electron 
and hole circles around the k=0 point. 
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